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ABSTRACT 

Primordial magnetic fields possibly generated in the very early universe are one of 
the candidates for the origin of magnetic fields observed in many galaxies and galaxy 
clusters. After recombination, the dissipation process of the primordial magnetic fields 
increases the baryon temperature. The Lorentz force acts on the residual ions and 
electrons to generate density fluctuations. These effects are imprinted on the cosmic 
microwave background (CMB) brightness temperature fluctuations produced by the 
neutral hydrogen 21cm line. We calculate the angular power spectrum of brightness 
\& , temperature fluctuations for the model with the primordial magnetic fields of a several 

nano Gauss strength and a power-law spectrum, ft is found that the overall amplitude 
C"| . and the shape of the brightness temperature fluctuations depend on the strength and 

the spectral index of the primordial magnetic fields. Therefore, it is expected that 
q , the observations of the CMB brightness temperature fluctuations give us a strong 

!_h ' constraint on the primordial magnetic fields. 

Key words: cosmology: theory - magnetic fields - large-scale structure of universe 



1 INTRODUCTION 

Observations reveal existence of magnetic fields on very large scale, i.e., galaxies and clusters of galaxies. It is found that these 
magnetic fields typically have a few uGauss strength s and relatively large coherent scales, i.e., a few tens of kpc for clusters of 
galaxies and a few kpc for galaxies iKronberslll994l) . The origin of such magnetic fields has not yet known while many ideas 
have been proposed. Perhaps it may be one of the most important remaining problems of cosmology and astrophysics to find 
out the origin and evolution of magnetic fields in the history of the universe. 

The most conv entional idea is that such magnetic fields were formed due to the astrophysical processes such as Biermann 
battery l|Biermannlll95(t in stars and supernova explosions. Then these seed magnetic fields were amplified by the dynamo 
process. Eventually, supernova winds or active galactic nuclei (AGN) jets may spread these magnetic fields into inter-galactic 
medium (for a comprehensive review see IWidrowlEooi) . However, it is still little known about the efficiency of the dynamo 
process in the expanding universe. It is particularl y difficult for the astro physical processes to explain observed magnetic fields 
with very large coherent scales in galaxy cluster s kim et al.lll99d.ll99llh A recent observation suggests existence of magnetic 
fields in high redshift galaxies ijKronberg et alJll992l) . These galaxies may be dynamically too young for the dynamo process 
to play a role. 

An alternative scenario is that magnetic fields were formed in the very early universe. Many authors have suggested 
various generation mechanisms of the primordial magnetic fields in the early universe. One can introduce an exotic coupling 
between electro-magnetic and scalar fields to generate magnetic fields during the inflation epoch. Or one can consider bubble 
collisions d uring cosmologica l phase transitions such as QCD or electroweak for generation of magnetic fields. For a detailed 
review, see lQiovanninil ((2004) . In this alternative scenario, one may directly obtain the nano Gauss primordial magnetic fields 
which are sufficient enough to explain /iGauss magnetic fields observed at present since the adiabatic compression due to the 
structure formation could easily amplify the primordial magnetic fields by a factor ~ 10 3 . In this case, there is no need of the 
dynamo process. However, if the seed magnetic fields generated in the early universe were too weak, the dynamo process is 
required even in this scenario while the coherent length could be very large unlike the astrophysical processes. 
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If the magnetic fields were generated in the very early universe, such primordial magnetic fields may play an important 
role for various cosmological phenomena. There are many previous works to constrain the strength of the primordial magnetic 
fields from Big Bang Nucleosynthesis (BBN), temperature anisotropies and polarization of Cosmic Microwave Background 
(CMB), or structure formation. These constraints give us clues to the origin of large-scale magnetic fields, and when and how 
magnetic fields were generated. 

The primordial magnetic fields affect on BBN through the enhancement of the cosmological expansion rate as neutrinos 
or the modification of the reaction rates of the light elements. The limit on the magnetic field st rength from BBN is Bo < 7 x 
10 Gauss where Bo is the comoving magnetic field strength Jcheng et alJll99d iKernan et allll996h . 

The black-body energy spectrum of CMB could be distorted by the primordial magnetic fields through the dissipation 
process. The severe constraint on the spectrum disto rtion from CO BE/FIRAS observation leads to Bo<3 x 10 _8 Gauss on 
comoving scales of 400-600pc at redshift z > 2 x 10 6 jjedamzik et al-lboooh . 

The primordial magnetic fields produce CMB temperature anisotropies on various scales. On very large scales, coherent 
magnetic fields generate the anisotr opic expansion of th e universe. Magnetic fields with the present-horizon-size coherent scale 
are constrained as Bo ^ 10 -9 Gauss ( B arrow et al Jl997l) . On intermediate or small scales, magnetic pressure modifies the acous- 
tic oscillations of baryon-photon fluids. Or Alfven mode (fluid vorticity) induces temperature anisotropies. Constraints on the 
primordial magnetic fields with lOOMpc-lMpc coherent scales from WMAP and other CMB ex periments are Bo ^ 10 -8 Gauss 
llAdams et Zlll996t iMack et al-lEooa LewislboO^ lYamazaki et al J I2OO5I: iTashiro et al.ll2005l) . Moreover CMB polarization 
suffers from magnetic fields. The Faraday rotat ion of the polarization is caused when CMB crosses the ionized inter-galactic 
medium (IGM) if there exist magnetic fields JScoccola et alJ I2OO4I : [pCosowskv et alJ I2OO5I) . The fluid vorticity induced by 
magnetic fields can generate B-mode (parity odd) polarization as well as E-mode (parity even) JSubramanian et al.l 120031 : 
ITashiro et aH Bo05). We expect to have much stringent limits on the primordial magnetic fields by future CMB observations, 
e.g., Planck (http://www.rssd.esa.int/Planck). 

The primordial magnetic fields may play an important role for the structure formation in the universe. The Lorentz 
force of t he primordial magnetic fields could induce density fluctuations once the universe becom es transparent after recom- 
bination iWassermanlll97a : iKim et alJll996l : Fsubramanian fc Barrowlll99gl: iGopal fc Sethill2003h . The magnetic tension and 
pressure are more effective on small scales where the entanglements of magnetic fields are larger. Therefore, if there exist the 
primordial magnetic fields, it is expected that there is the additional power in the density power spectrum o n small scales 
and these power induces the early structure formation JSethi fc Subramanianll2005t ITashiro fc Sugivamall200r?) . 

Thermal evolution of baryon s after decoupling from pho tons, redshift z ~ 200, could be also modified by the existence 
of the primordial magnetic fields JSethi fc Subramanianl Bo05l since the dissipation of the primordial magnetic fields work as 
the heat source. The dissipation mechanisms are the ambipolar diffusion and the direct cascade decay of magnetic fields. 

Both the structure formation and the thermal evolution of baryons are closely connected to reionization of IGM. Hence, 
we can suspect that the reionization process can be strongly affected by the existence of the primordial magnetic fields. 

One of the best probes of the reionization process of IGM is the CMB brightness temperature fluctuations induced by neu- 
tral hydrogens through the 21cm line. In this paper, therefore, we examine the effect of the primordial magnetic fields on the 
CMB brightness temperature fluctuations produced by the redshifted hydrogen 21cm line. The hydrogen 21cm line is caused 
by a spin flip of the electron in a neutral hydrogen atom. The neutral hydrogen atoms produce the brightness temperature 
fluctuation of CMB by 21cm line absorption from and emission into CMB. The amplitude of the brightness temperature fluc- 
tuations depends on the hydrogen density, ionization fraction of hydrogens and hydrogen temperature. Therefore observations 
of the brightness temperature fluctuations at wave length 21(1 + z)cm reveal the density fluctuat ions and the ionization pro- 
cess at redshift z jLoeb fc ZaldarriagalEool iBharadwai fc Alilbooi ICooravlEool lAli et aljEool) . Recent efforts to measure 
redshifted 21cm line such as LOFAR (http://www.lofar.orgl, MWA (http://web.haystack.mit.edu/arrays/MWA/MWA.htmll 
and SKA (http://www.skatelescope.orgl, will soon give us clues of dark ages of the universe. 

This paper is organized as follows. In Sec. II, we discuss the effect of the primordial magnetic fields on the hydrogen 
temperature and the density fluctuations. In Sec. Ill, we summarize CMB bright temperature fluctuations produced by the 
hydrogen 21cm line. In Sec. IV, we compute the angular power spectrum of bright temperature fluctuations with the primordial 
magnetic fields and discuss the effect of the primordial magnetic fields. Sec. V is devoted to summary. Throughout the paper, 
we take WMAP val ues for the cosmolog ical parameters, i.e., h = 0.71 (H = hx lOOKm/s ■ Mpc), T = 2.725K, ft 2 f2 b = 0.0224 
and /i 2 ft m = 0.135 iSpergel et alj F2003). And Ti and c are Planck's constant over 2-7T and speed of light, respectively. 



2 EFFECTS OF THE PRIMORDIAL MAGNETIC FIELDS ON THERMAL HISTORY AND 
STRUCTURE FORMATION 

After recombination, the primordial magnetic fields affect the thermal history and the structure formation of the universe. In 
this section we summarize how the existence of magnetic fields modifies them. 

Let us discuss the evolution of the primordial magnetic fields and the resultant power spectrum. First, we postulate that 
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the primordial magnetic field lines are frozen-in to baryon fluid in the early universe. This assumption leads to the time 
evolution as 

B (x) =a 2 (t)B(t,x), (1) 

where a(t) is the scale factor, which is normalized to the present value, and Bo is the comoving strength of magnetic fields. 
Since the baryon fluid is highly conductive through the entire history of the universe, this relation can be hold as far as 
magnetic fields do not suffer from non-linear processes such as direct cascade induced by turbulent motion of eddies. 

Next, we assume that the primordial magnetic fields are isotropic and homogeneous Gaussian random fields, and have a 
power-law spectrum as 

(B (M (ki)B,y(k 2 )) = ^<5(k! - k 2 ) - Bl{k), (2) 

4nk 3 B%(k) = (n + 3) (±) B 2 for (k < k c ), (3) 

where k c is the cutoff wave-number and n is the spectral index. Here B c is the root-mean-square amplitude of the magnetic 
field strength in real space. 

The cutoff in the power spectrum appears due to dissipation of magnetic fields associate d with the nonlinear direct 
cascade process Jjedamzik et alJll998tlSubramanian fc Barrowlll'998l : lBaneriee fc Jedamzikll2004F) . The direct cascade process 
transports the magnetic field energy from large scales to small scales through the breaking of flow eddies and generates the 
peak in the energy spectrum resultantly. The time-scale of the eddy breaking at the scale I is l/v, where v is the baryon 
fluid velocity. The direct cascade occurs when the eddy breaking time-scale is equal to the Hubble time _ff _1 . Once the direct 
cascade takes place, the "red" power-law tail is produced in the power spectrum of magnetic fields. In order to simplify 
following calculations, however, we assume that the direct cascade process produces a sharp cutoff instead of the power-law 
cutoff in the spectrum. 

Once the universe becomes transparent after recombination, baryons decouple from photons and their velocity starts to 
increase. Eventually the velocity achieves the value determined by the equipartition between the magnetic field energy and 
the kinetic energy of the baryon fluid, namely the Alfven velocity, va = cBo / \J 'in pboa(t) where pb is the baryon density and 
the subscript denotes the present value. Accordingly, the comoving cutoff scale induced by the direct cascade process is 
given by 

Ha 52MPC- f^^r 1 (^- V/2 X : ' 



kc « 2tt— « 52MPC" 1 c -^2- . (4) 

va VInGauss/ ^ 0.0224 J ^0.135/ w 

Note that this comoving cutoff scale is time independent in the matter dominated epoch. 



2.1 Dissipation of magnetic fields 

The energy of the primordial magnetic fields is dissipated by the ambipolar diffusion and the direct cascade process. The 
dissipation of magnetic field energy gives a considerable effect to the thermal evolution of hydrogens. If the energy of the 
primordial magnetic fields with the strength Bo is instantaneously converted to the thermal energy of hydrogens at redshift 
2, the hydrogen temperature Tk can be estimated as Tk ~ 30(1 + z)(£>o/lnGauss) 2 K. 

The ambipolar diffusion is caused by the velocity difference between ionized and neutral particles. Because magnetic 
fields accelerate ionized particles by the Lorentz force while they give no effect on neutral particles, the difference of velocity 
between ionized and neutral particles arises. This difference induces the viscosity of ionized-neutral baryon fluid. Accordingly, 
the energy of magnetic fields is dissipated into the fluid. 

On the other hand, the direct cascade is the nonlinear process which induces coupling between different modes. When 
the time-scale of the eddy breaking equals the Hubble time, the magnetic field energy is transported from large scales to 
small scales by the nonlinear process and the transported energy ends up dissipated. The dissipation process due to the direct 
cascade, therefore, shifts the cutoff scale of the spectrum of the magnetic fields to larger scales. However it is shown in Eq. |g) 
that the cutoff scale does not evolve during the matter dominate d epoch. So we can conclude that the direct cascade is not 
the major source of the dissipation of the magnetic fields. In fact. lSethi fc Subramanianl ^OOflh have investigated the effect of 
these dissipations on the cosmological thermal history and have showed that the effect of the ambipolar diffusion dominates 
that of the direct cascade process. In this paper, therefore, we take into account only the ambipolar diffusion as the dissipation 
process. 

The energy dissipation rate of the ambipolar diffusion is described as JCowlinelll95£l) 

r = i6^: 1(VxB)xB|2 ' (5) 
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where x e is the ionization fraction and \ i s the drag coefficient for which we adopt the value computed bv lDrairie" et al l dl983l) 
: x — 3.5 x 10 13 cm 3 g _1 s _1 . For the Gaussian statistics of the primordial magnetic fields, Eq. the energy dissipation rate 
can be rewritten as 

r = 192^^1° / dkl J dk,Bl{k,)Bl(k 2 )klki (6) 
The evolution of the hydrogen temperature with the ambipolar diffusion, which is described by F, is given by 

# = + T^T 2 ^^ - ^ + ttt — > (7) 

at a 1 + x e im e c Lok^n^ 

where p 7 , <tt, m e , T 7 and A?b are the CMB energy density, the Thomson cross-section, the electron mass, the CMB temperature 
and the Boltzmann constant, respectively. The dot represents a derivative with respect to time. For the standard cold dark 
matter dominated universe model, the residual ionization fraction after recombination is x £ « 10" 4 until the universe rciomzcs. 
However, it is possible that higher hydrogen temperature due to the magnetic field dissipation makes the collisional io nization 
effect ive and resultantly, the ionization rate becomes higher. The evolution of the ionization fraction is described as ((Peeblesl 



^j- = [/3 e (l - x e ) - a e n h x 2 e ] C + 7 e n b (l - x e )x e , (8) 
where 

_ m e fc B T k 3/2 / AE \ 

is the ionization rate out of the ground state with the ground state binding energy AE = 13.6eV, a e is the recombination 
rate to excited states and C is a suppression factor (for details , see |Peebleslll96Sl ll993h . The 7 e in the second term is the 
collisional ionization rate and we utilize the fitting formula bv IVororlo*Xi[l99'it) . The collisional ionization is suppressed by 
exp(— AE/knly). Therefore, the collisional ionization does not dominate the first term until the hydrogen temperature exceeds 
10 5 K. 

We calculate Eqs. J7J and JHJ by using the modified RECFAST code JSeaeer et alJll999h . We take into account the 
ambipolar diffusion as the dissipation process of magnetic fields and ignore other heat sources, e.g., heat from stars and 
galaxies. We plot the hydrogen thermal evolution in Fig. Q In the redshift higher than z ~ 200, the hydrogen temperature 
traces the CMB temperature due to the second term in the right-hand side of Eq. J7J , which represents the Compton scattering. 
After z ~ 200 hydrogens have decoupled from photons. Since then, the hydrogen temperature declines faster than the CMB 
temperature as Tk oc (1 + z) 2 . The heating by the dissipation of the magnetic field energy gradually becomes effective and 
eventually the hydrogen temperature goes up. The heating efficiency of the dissipation depends on not only the magnetic 
field strength but also the power-law index of the spectrum. This is because the magnetic fields with a shallower spectrum 
produce the velocity difference between ionized and neutral hydrogens on a wider range of scales. Therefore, the dissipation 
from magnetic fields with a smaller spectral index is more effective than that with a larger spectral index since we normalize 
the magnetic strength at the cutoff scale. 

Fig-IUshows the evolution o f the ionization fraction. If t he hydrogen temperature is higher than 5 x 10 5 K, the collisional 
ionization becomes effective fsee lSethi fc Subramanianll2005l) . In our case, however, the hydrogen temperature never exceeds 
this value so that the modification of the ionization fraction by the magnetic field dissipation is very little. 



2.2 Generation of density fluctuations 



Primordial magnetic fields generate the density fluctuat ions after recombination iWassermanll978tlKim et al.ll99rllSubramanian fc Barrowl 
19981: iGopal fc Sethilbood ISethi fc SubramanianlEooi) . The magnetic tension and pressure are more effective on small scales 



because the entanglements of magnetic fields are larger. Therefore, on small scales, the additional density power by magnetic 
fields is expected to dominate the primordial density power spectrum produced by inflation. The evolution equations of the 
density fluctuations with the primordial magnetic fields are described as, 

^ = ~ 2 llft + 47TG( -P hSh + Pdm^dm) + S(t, X), (10) 

s(t, x) ^ v -« v> : Bo w^ Bo W), (id 

at- 1 a at 

where pdm is the dark matter density, and 5b and <5dm are the density contrasts of baryons and dark matters, respectively. 
In order to solve these equations, we define the total matter density p m and the matter density contrast <5 m as 
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Figure 1. The evolutions of the hydrogen temperature. The dashed, dashed-dotted and solid lines are the hydrogen temperature for 
magnetic fields with B c = 3nGauss and n = 0, B c = InGauss and n = —2 and B c = InGauss and n = 0. The dotted line is the hydrogen 
temperature in the standard cosmology (without magnetic fields) which is proportional to (1 + z)~' 2 once hydrogens decoupled with 
photons at z ~ 200. The thin solid line represents the CMB temperature which is proportional to (1 + z) • 
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Figure 2. The evolution of the ionization fraction. The dashed, dashed-dotted and solid lines are the ionization fraction for magnetic 
fields with B c = 3nGauss and n = 0, B c = InGauss and n = —2 and B c = InGauss and n = 0. The dotted line is the ionization fraction 
in the standard cosmology (without magnetic fields). In the neutral hydrogen number density, there is no significant change due to the 
energy diffusion of magnetic fields. 




p m = p b + pdm, (13) 

_ (Pb<5b + Pdm<5dm) 

Om = • 

Pm 

From Eqs. 1101 and (11 2^ . the evolution of <5 m is given by 



(14) 



„, = -2-—— + 47rGp m 5 m H S(t, x). (15) 

at 1 a at p m 

The solution of Eq. JT^J can be acquired by the Green's function method, 

8 m = AWD 1 (t)+B(x)D 2 (t)-—D 1 (t) dt w{tl) +^D 2 (t) ^ dt ^ , (16) 

where D\ (t) and D2 (t) are the homogeneous solutions of Eq. 1151 and W is the Wronskian and is expressed as 

W(t) = Dx{t)Di(t) - Lh{t)Di(t), (17) 

and ti denotes the initial time. 

The first and second terms of Eq. I|16fl correspond to the growing and the decaying mode solutions of primordial density 
fluctuations and the third and fourth terms are the ones generated by the primordial magnetic fields. We represent the former 
two terms as <5 m p and the latter twos as <5 m M. Here we only consid er th e growing s o lution for 5 m p. The analytic solution of <5 m M 
in the matter dominated epoch is obtained bv I Wassermanl lll97ct) and lKim et alJ dl996h . In this paper, we concentrate on the 
evolution of fluctuations before reionization z > 10 when the universe is matter dominant and the contribution from the dark 
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energy is negligible. Therefore we can utilize their solution. In the matter dominated epoch, Di(t) oc t 2 / 3 and D2(t) oc t 
Accordingly the terms generated by magnetic fields of Eq. 1 Kit can be written as 



SmM = 



(5) 



2/3 



ttS(ti,x). 



(18) 



Once the density fluctuations were generated by the primordial magnetic fields, they grow due to the gravitational instability 
so that the growth rate is same as the primordial density fluctuations. 

Next, we calculate the power spectrum of the matter density fluctuations. Assuming that there is no correlations between 
the magnetic fields and the primordial density fluctuations for simplicity, we can describe the matter power spectrum as 

P m (k) = P mP (fc) + P mM (fc) = <|<W(fc)| 2 ) + (|<W(*0| 2 ), (19) 

where S nl p(k) and <5 m M(fc) are the Fourier components of <5 m p and <? m M, respectively, and ( ) deno tes the ensemble average. 

We numerically calculate P m p(fc) by using the CMBFAST code iSeliak fc Zaldarriagalll99a) . while P m M(fc) is obtained 
from Eq. (ITS! as 

,n 2 



-PmM(fc) = 

where 



(fim) V 



if 



V 4:np b0 a 3 (ti) 



(t\ 2 ''i 3/t\- 1 



7 2 (fc), 



I 2 {k) = <|V ■ (V x Bo(x)) x B (x)| 2 ). 

Applying isotropic Gaussian statistics Eq we can rewrite the nonlinear convolution Eq. 112 It as 



I 2 (k) = 



J dkl J ^ ^(fcOPgCk-k.l) [2fc5fc 3 M + fc 4 fc 4 (1 _ ^ 



■2k? kin 3 ] 



(20) 



(21) 



(22) 



where fi is /J, = k • ki/fc&i. The integration of I 2 (k) is determined by the value of integrand at fei = k c and |ki — k| = k c 
because the power spectrum B^(k) has the power law shape with sharp cutoff at k — k c . Note that the direct cascade process 
of the magnetic fields from large scales to small scales likely produces a power-law tail instead of sharp cutoff in the power 
spectrum above k c as we discussed before. However, it is still true that most of the contribution on the integration I 2 (k) 
comes from the peak of the spectrum at k c as well as the sharp cutoff case, although some corrections may be needed. 

Let us introduce an important scale for the evolution of density fluctuations, i.e., the magnetic Jeans length. Below this 
scale, the magnetic pressure gradients, which we do not take into account in Eq. (I15L counteract the gr avitational force and 
prevent further evolution of density fluctuations. The magnetic Jeans scale kMj ~ v\/H is evaluated as iKim et alJll996h 

R v -1 / tJn \ /l2 ^ \ X / 2 



KMJ = 07T = 12. 



So 



7Mpc _1 ( — 
\ In 



Gauss / 



ft 2 fi b 
0.0224 



0.135 



(23) 



Here, the primordial magnetic fields have a power-low spectrum with spectral index n and sharp cutoff at k c so that the 
magnetic Jeans scale can be written from Eqs. |gj and i'l'At as 



(24) 




&MJ = Ti — \ 7T 



Baryon density fluctuations below the magnetic Jeans scale (the mode k > kuj) oscillate and do not grow. 

Using Eq. 122H . we numerically calculate the matter power spectrum P m (fc). We show the evolution of fc 3 P m (fc) for the 
primordial magnetic fields with B c — InGauss and n — 1 in Fig. El The contribution from the density fluctuations generated 
by the primordial magnetic fields is dominated on small scales ( < lOOkpc). Since the growth rates of both P m p and P m M are 
t i / 3 oc 1/(1 + z) 2 , the amplitude of the total matter power spectrum P m is proportional to 1/(1 + z) 2 and the comoving scale 
on which P m M starts to dominate P m p stays constant (k ~ 15Mpc _1 ) as is shown in Fig. 01 

In Fig. 21 we plot the matter power spectra for models with different magnetic field amplitudes at z = 30. We can 



in the limit of k/k c <C 1 a s I (k) ~ aB, 



2n + 10£.2n + 7 



+ [3Bj.k 4 where a and f3 are coefficients 



analytically estimat e Eq. 

which depend on n jKim et al.lll99fil : lOopal fc SethlbOQ.^ . Here we employ the fact that the cutoff scale k c is proportional 
to P^T 1 as is shown in Eq. ||1J. The former term dominates if n < —1.5, while the latter one dominates for n > —1.5. 
Accordingly, the matter power spectrum, k 3 P m (k), is proportional to B 2n+10 k 2n+7 for n < —1.5 or to Pjfc 4 for n > —1.5. 
These dependences of k 3 P m (k) on B c and n can be found in Fig. 0] 



3 BRIGHTNESS TEMPERATURE FLUCTUATIONS BY THE 21CM LINE 

In this section, we discuss the calculation of the CMB brightness temperature fluctuations generated by the hydrogen 21cm 
line and the angler power spectrum of them. First, let us discuss the spin temperature on which the amplitude of the brightness 
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Figure 3. The redshift evolution of the matter power spectrum for the model with B c = InGauss and n = 0. The dashed, dotted and 
solid lines represent the matter power spectrum at z = 20, z = 30 and z = 40 (from top to bottom). The thin solid line is the matter 
power spectrum without magnetic fields at z = 40. 
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Figure 4. The matter power spectra for models with different magnetic field strengths. The dashed, dotted and solid lines are the power 
spectra with the primordial magnetic fields B c = 3nGauss and n = 0, B c = InGauss and n = —2 and B c = InGauss and n = (from 
top to bottom) at z = 30. 



temperature fluctuations depends. The spin temperature is defined through the ratio of the number density of hydrogen atoms 
between in the excited state and in the ground state of the 21cm transition, 

^=3exp(-^), (25) 
no V Is J 

where subscripts 1 and denote the excited and ground states and T» corresponds to the energy difference between these 
states, T* = 0.0682K. The factor 3 in Eq. 

comes from the ratio of the spin degeneracy factors. The evolution of the 
number density of hydrogen atoms in the ground state is described as 

^ - r-7 — no = - jz , \„, , [-n (C*oi + B 01 I V ) + m (C 10 + A 10 + B 10 J„)] , (26) 
Oz l + z (l + z)H(z) 

where Coi and Cio are the collisional excitation and de-excitation rate coefficients, Aio is the Einstein A-coefficient, Boi and 
Bio are the Einstein B-coefRcients and I v is the specific intensity of CMB at frequency v. The collisional de-excitation coefficient 
is written as Cio = 4/tio(Tk)nH/3 where nu = np + ni is the neutra l hydrogen number density and Kio(Tk) is tabulated as a 
function of the hydrogen temperature Tk jAllison fc Dalgarnolll969l) . Typically Kio(Tk) is order 10 10 cm 3 s 1 in Ty, < 1000K. 
The excitation coefficient is related with the de-excitation coefficient by the detailed balance as Cqi = C\o exp(—2nhv / kTy) . 
The Einstein coefficients are coupled by the Einstein relations as 



B 01 = 3Bi =^ i ^Mio, r27) 



where Aio = 2.85 x 10 15 s . The evolution of the number density of the excited hydrogen atoms is also written as 
dni 3 

Th^TTz'"'^ (l + z)H{z) 

From Eqs. JHJ and we can derive the evolution of the spin temperature 



~ ni = ~TT, — \WT\ I n ° ( Co1 + Boi^) - ni (Cio + A 10 + Bi I v )] . (28) 
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d_ 

dz 



(l + z)H(z) 



C10 + 



T, ' 



(29) 



where we employ the approximation exp(— T»/T) ~ 1 — T*/T, and the Rayleigh- Jeans law, /„» = 2vlkB.T 1 / (? with f» = 
fesT* / '(2irh) because we are interested in the epoch when T s , Tk, T 7 2> T». 

The emission and absorption of the 21cm line from neutral hydrogens affect the CMB brightness temperature at the 
frequency v*/{l + z) by 



Tk 



T 



+ / I- + : 

where n is a unit direction vector and r is the optical depth of the 21cm line, 
3nc 3 nnhAio 1 



ST b 



16u'ik B T s H{z)' 

The fluctuations of the brightness temperature at the frequency v*/(l + z) are obtained from Eq. %MM as 

8% 



(1 



1 



(T s — T 7 ) 5 m + T~i 



(30) 



(31) 



(32) 



where we ignore the fluctuations of the CMB temperature and the ionization fraction and we assume that the hydrogen number 
density contrast is corresponding to the matter density contrast. From Eq. 1321 . it is found that the brightness temperature 
fluctuations at the frequency + z) are determined by the density and the spin temperature fluctuations at redshift z. 

Note that we ignore the effect of the line-of-sight component of the neutral hydrogen's peculiar velocity. The line-of-sight 
component of the peculiar velocity affects the optical depth through the distortion of th e redshift space. The c ontribution 
from the peculiar velocity could be substantially large on very large scales (for details, see iBharadwai fc AlilbooJ . while it is 
negligible on scales comparable to the future observations' survey area. Therefore we ignore the term of the peculiar velocity. 

In order to calculate Eq. 13211 . we need to know the time evolution of the spin temperature fluctuations. The evolution 
equation can be obtained from Eq 12911 as 



d_ 

dz 



5% 
T 



4T 



(Cva , A W \5T S fl 1 \ <5n H // 1 1 



dlnn 1 \ 6T k 
dlnT k T k J 10 T k 



H(z)(l + z) LV T k ' T« J T B ' VT k T s y^ ±u n H ' VVT k TjdlnT^ T^J " iu T k J ' (33) 

where we ignore the fluctuations of the CMB temperature and the ionization fraction. Because of the energy pumping from 
CMB, strictly speak ing, the ordinary adiabat ic relation, i.e. 5T k /T k = (7— l)5nu/nu where 7 is adiabatic index and 7 = 5/3, 
is broken. However, IBharadwai fc Alii J2004) have shown that the adiabatic relation is still valid until the star formation 
process proceeds. Hence, hereafter, we assume <5T k /T k = (7 - l)5n H /n H in Eq. 

Let us calculate the angular power spectrum of brightness temperature fluctuations Ci(u*/(1 + z)). We can expand 
the brightness temperature fluctuations 5T b {v*/{l + z), n) into spherical harmonics with the expansion coefficients ai m . 
Accordingly the angular power spectrum C;(i/*/(l + z)) = ^|a; m | 2 ^ can be written as 



C, 



(1 + *) 



\air, 



2 ) =4tt 



dk* 
(2tt) 3 



|<5T b 



(1 + ^) 



k 



3i (Hvo-V(z))), 



(34) 



where ji(x) is the spherical bessel function, r\ is the conformal time and &T^(v,k) is the Fourier component of <5Tb at the 
frequency v which is calculated by Eq. I32H . 



4 RESULTS AND DISCUSSION 

We can calculate the angular power spectrum of the brightness temperature fluctuations produced by the 21cm line with the 
existence of the primordial magnetic fields from Eq. 1341 . In this section, we show the evolution of the spin temperature and 
the angular power spectrum for three different primordial magnetic fields, i.e., (B c , n) — (InGauss, 0), (InGauss, — 2), and 
(3nGauss, 0). 

First we discuss the effect of the primordial magnetic fields on the spin temperature. Fig. Q3shows the redshift evolution 
of the spin temperature. The spin temperature is well approximated to the steady state of Eq. I|29^ as l)Fieldlll959h 

T S = Z±&, (35) 
1 + y 

where y is the collisional efficiency and written as 
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The behaviors of the spin temperature in Fig. can be explained by Eq. 1351 . In the high redshift (z > 100), because the 
collisional efficiency Cio is high (y 2> 1) so that the spin temperature couples the hydrogen temperature, i.e., T s ~ T k - As the 
universe expands the collisional efficiency becomes low (y <C 1) so that the spin temperature approaches the CMB temperature 
T s ~ T 7 . The primordial magnetic fields make the hydrogen temperature increase due to the dissipation. The collisional 
efficiency increases as the hydrogen temperature increases. Accordingly the combination j/Tk in the numerator of Eq. I|350 
increases although y decreases. Therefore the spin temperature becomes higher than the one without the magnetic fields. If 
the primordial magnetic fields are strong enough, the hydrogen and spin temperatures can exceed the CMB temperature as 
is shown in Fig- El 

Next, let us show the angler power spectra of the brightness temperature fluctuations by the 21cm line for the various 
primordial magnetic fields in Fig . El For reference , we p lot the angler power spectra without magnetic fields, which are 
consistent with previous works bv lBharadwai fc Alii J2004) . 

The primordial magnetic fields affect the overall amplitude of the angular power spectrum and the shape of the angular 
power spectrum on large Vs. The modification on the shape of the angular power spectrum is caused by the density fluctuations 
generated by the magnetic fields on small scales. From Eq. 13211 . the brightness temperature fluctuations are determined by 
the density fluctuations <5 m and the spin temperature fluctuations ST S /T S . Applying the steady state approximation and j<1 
to Eq. 1)33^. we can obtain 8T S /T S in the leading order of y as 



5% = (T\_ 
T s V \T-, 

where we ignore the temperature derivative term of k and substitute T s with Eq. 1351 . Accordingly Eq. 13211 leads to 



(37) 



8T h = -—y [2T k - (1 + 7 )T 7 ] <5 m . (38) 
1 + z 

It turns out that the brightness temperature fluctuations are simply proportional to the density fluctuations. Therefore the 
additional "blue" density power spectrum produced by the primordial magnetic fields on small scales appears as the "blue" 
angler power spectrum on large I's. The magnetic Jeans scale corresponds to £mj ~3x 10 5 (lnGauss/_B c ). The angular power 
spectrum C; increases until Imj and starts to show the oscillatory behavior above I > Imj- 

The evolution of the power spectrum amplitude with the primordial magnetic fields is classified into two regimes: T 7 > T s 
and T-y <T S . In T 7 > T s , neutral hydrogens absorb the CMB photons while neutral hydrogens emit the 21cm line into CMB 
in T 7 <T S . If there is no magnetic fields, T 7 > T s until the formation process of stars and galaxies takes place. On the other 
hand, the dissipation of the primordial magnetic fields makes the hydrogen temperature higher than the CMB temperature 
and T B can exceed T 7 at high redshift as mentioned above. 

In the regime of T 7 > T a , the amplitude of the angular power spectrum with the primordial magnetic fields declines as is 
the case with no magnetic fields due to the fact that T s approaches to T 7 . However, the decline of the power spectrum is more 
prominent for the model with the primordial magnetic fields (see bold and thin solid lines on the top and middle panels of 
FIG.EJ- The reason is following. For the model without magnetic fields, T 7 ^> Tk. From Eq. I|38|l . therefore, the angular power 
spectrum Ci oc |(1 +7)T 7 | 2 . If the primordial magnetic fields exist, on the other hand, Tk increases due to the dissipation. 
Accordingly Ci oc |(1 + 7 )T 7 - 2T k | 2 < |(1 + t)T 7 | 2 as far as T k < (1 + Therefore the suppression of the angular power 

spectrum of the model with the primordial magnetic fields is more prominent. 

Once the hydrogen temperature T k exceeds the CMB temperature T 7 (or equivalently T s > T 7 ), the amplitude of 
the brightness temperature fluctuations starts growing as the universe evolves. In the low redshift universe, the hydrogen 
temperature becomes high enough that we can assume the CMB temperature is negligible to the hydrogen temperature, 
T 7 <C Tk. Hence the fluctuations of the brightness temperature can be approximated as 

ST h = -J—2yT k S m . (39) 
1 + z 

Perhaps one might think that the amplitude of the magnetic fields can be determined by measuring the angular power 
spectrum since Tk depends on the amplitude. However it is not the case. Because y is related to the hydrogen temperature 
through Cio/Tk, the overall dependence of 5Tb on the hydrogen temperature comes from Cio- In 8T^ > 100K, Cio is not a 
strong function of the hydrogen temperature. Accordingly, the amplitude of the angular spectrum does not depend much on 
the strength of the primordial magnetic fields as is shown in Fig. El Therefore, in the regime of T 7 < T s , we can only measure 
the strength of the primordial magnetic fields through the turnover of the spectrum on high Z's. 



5 SUMMARY 

In this paper we investigate the effect of the primordial magnetic fields on the CMB brightness temperature fluctuations 
produced by the hydrogen 21cm line. 

The brightness temperature fluctuations depend on the density fluctuations, the hydrogen temperature and the ionization 



10 Tashiro, H. et al. 




Figure 5. The evolution of the spin temperature. The solid lines, the dotted and dashed lines denote the spin temperature, the hydrogen 
temperature and the CMB temperature, respectively. We show these temperature for three different magnetic field spectra: B c = InGauss 
and n = 0, B c = InGauss and n = —2 and B c = 3nGauss and n = (from top panel to bottom panel). For reference, we plot the spin 
temperature and the hydrogen temperature in the standard cosmology with no primordial magnetic fields as the thin solid lines and 
dotted lines. In all cases, the spin temperature closes in the hydrogen temperature in high redshift and approaches the CMB temperature 
asymptotically. 

fraction whose evolutions are affected by the primordial magnetic fields: the dissipation of the magnetic field energy works 
as heat source and the Lorentz force generates the density fluctuations after recombination. The primordial magnetic fields 
with nano Gauss strength can heat hydrogens up to several thousands degrees Kelvin by the ambipolar diffusion and produce 
dominant additional "blue power" in the matter density spectrum on scales smaller than lOOkpc. 

Through these effects, the primordial magnetic fields modify the angular power spectrum of the brightness temperature 
fluctuations. First, an additional blue spectrum can be found due to the existence of the blue power in the matter density 
spectrum. Secondly, the overall amplitude of the angular spectrum can be modified by the magnetic fields. If the hydrogen 
temperature is lower than the CMB temperature, the amplitude becomes smaller than the one without magnetic fields. The 
difference depends on the hydrogen temperature. On the other hand, if the hydrogen temperature becomes higher than the 
CMB temperature, the amplitude of the angular power spectrum, which depends little on the magnetic field amplitude in 
this case, exceeds the one without magnetic fields. 

Perhaps each feature of the brightness temperature fluctuations is not a direct evidence for the existence of the primordial 
magnetic fields. Possible existence of the isocurvature mode on small scales can also boost the spectrum on small scales, while 
the higher hydrogen temperature might be achieved by the other heat source, e.g., the decaying dark matter. However, both 
blue spectrum and modification of overall amplitude of the angular power spectrum may provide a unique evidence for the 
existence of the primordial magnetic fields. Or at least we can set a stringent constraint on the shape and the amplitude of 
the primordial magnetic fields. 

In this paper, we ignore the reionization process due to the ordinary astronomical objects, i.e., stars, which provides 
a significant impact on the brightness temperature fluctuations, s ince we are focus o n the epoch before reionization. The 
reionization can be as late as z ~ 9 from the latest WMAP results JSpereel et alJl2006^ . Once the reionization process starts, 
the ionizing UV photons produced by stars make the hydrogen temperature high and the Ly-a pumping of the hydrogen 21cm 
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Figure 6. The angler power spectra of CMB brightness temperature fluctuations by the 21cm line at the frequency u*/(l + z) for 
z = 40, 30, 20, and 10 from the top to bottom panels, respectively. The dashed, dotted and solid lines represent the spectra for the models 
with B c = 3nGauss and n = 0, B c = InGauss and n = —2 and B c = InGauss and n = 0, respectively. The thin solid line in each panel 
is the spectrum of the model without the magnetic fields. 



transitions efficient. Accordingly, the coupling between the spin temperature and the hydrogen temperature becomes stronger. 
This reionization process could be also affe cted by the primordial mag netic fields due to the modification of the matter power 
spectrum which induces the star formation jTashiro fc Sugivamal2006l) . Therefore we need to consider the reionization process 
with the primordial magnetic fields in order to investigate the role of the magnetic fields below z < 10, which is beyond the 
scope of this paper. 
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